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Introduction 


Waves  in  a  viscoelastic  rod  have  been  studied  by  Morrison 
and  others  [1]  using  Laplace  transform  techniques.  The  form 
of  the  results  obtained  is  so  complex  that  it  is  difficult 
to  obtain  quantitative  information  without  extensive  numeri¬ 
cal  calculations;  nor  is  it  convenient  to  use  this  form  in  a 
Duhamel  type  integral.  For  this  reason  Bleich  and  Sackman 
[2,3]  have  proposed  an  exponential  interpolation  as  an  approx¬ 
imation  in  viscoelastic  wave  propagation  problems.  The  above 
papers  were  limited  to  the  consideration  of  materials  having 
only  one  typical  time  constant  like  the  standard  solid.  It 
is  the  purpose  of  this  note  to  extend  the  idea  of  this  expo¬ 
nential  interpolation  to  materials  which  have  a  more  complex 
response,  containing  two  or  more  time  constants.  As  before, 
[2],  the  interpolation  makes  use  of  the  known  discontinuity 
at  the  wave  front  and  the  known  long  time  value  of  the  re¬ 
sponse;  further,  it  is  anticipated  that  in  the  case  of  a  mate¬ 
rial  whose  elements  in  creep  or  relaxation  display  j  typi¬ 
cal  time  constants,  the  response  F(t)  in  a  wave  propagation 
problem  will  have  the  character  of  Fig.  1  and  may  be  approxi¬ 
mated  by  an  expression  of  the  form 

-t/t  -t/t  -t/t . 

F(t)  ~  F(«>)  -Ae  1  -he  2  -  , . .  -A.e  3  (l) 

where  t  is  the  time  after  arrival  of  the  first  signal  at 
the  location  considered. 

If  j  =  1  the  two  parameters  A  and  t  can  be  deter- 
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mined  from  the  values  P(0)  and  F(0)  which  are  easily  ex¬ 
tracted  from  the  Laplace  transform  solution  by  asymptotic 
methods,  but  this  method  is  insufficient  if  j  )>  2 .  The  pres¬ 
ent  note  approaches  the  latter  case  by  again  starting  with 
the  Laplace  transform.  The  transform  of  the  response  is  then 
expanded  in  a  series4,  the  inversion  of  which  gives  a  series 
consisting  of  j  exponential  terms,  c  J,  each  multi¬ 
plied  by  a  power  series  in  t.  Comparing  the  leading  terms 
of  this  expansion  with  Eq.  (1),  the  values  of  t ^  and 
are  obtained.  The  details  of  these  steps  are  illustrated  with 
the  following  example . 

Example 

Consider  the  problem  of  the  propagation  of  waves  in  a 
semi-infinite  (x  <  0)  rod  of  a  five-parameter  viscoelastic 
material,  initially  at  rest,  subjected  to  a  unit  stress  im¬ 
pulse  at  t  =  0,  x  =  0.  For  an  elastic  rod  with  Young's  mod¬ 
ulus  E  and  mass  density  p  the  stress  is  known  to  be 

c  =  H[t-x(p/E)^]  (2) 

where  H(t)  is  the  Heavy side  step  function.  Using  the  elas¬ 
tic-viscoelastic  analogy  [5]  the  Laplace  transform  of  the 
stress  in  a  viscoelastic  rod  subjected  to  similar  conditions 
is 

o  =  s"1  exp[-sx(p/E)?]  (3) 


* 


The  validity  of  the  expansion  has  been  discussed  in  [4], 


where  s  is  the  transform  parameter  and  E  is  a  rational 
function  of  the  transform  parameter.  For  the  five-parameter 
model  shown  in  Fig.  2 


s2(g  +G  +G  )  +  s[t-1(g  +G  )  +  t_1(G  +G  ) 1  +  t"1  t-1  G 

g  _  12  3 _ 1  2  3 _ 2  13  12  3 

(s+T-1)  (s+T-1) 

1  2 

(4) 

„  e  ( s+a )  (s+b) 

(s+c)  (s+d) 

where  a,  b,  c,  d,  and  e  are  known  constants.  The  inversion 
to  be  performed  is 

o  =  r>[s-eXp[-s*(£f*f$i}] 

Using  the  shift  theorem 

5  =  s'1  exp[-sx{p/e)i]  exp  [-sx(£  |±|  |±S)i  +  sx(p/e)5]  (5) 

o=H[t-x(p/e)i]  F  [t-x(p/e)i]  (6) 

where  F  is  well  behaved  as  s  ->  oo  ,  This  is  a  condition 
for  the  term  by  term  inversion  which  follows.  Using  the  ex¬ 
pansions 

00 

eX=n2oXn/nl  |x|<“ 


(1+x)2 


”  (-x)nr(n-l/2) 

n=0  r(-l/2)r(n+l) 


=  1  +  x/2  -x2/8  +  x3/l6  - 


I  x  |  <  1 


one  finds 

F(t)  =  £-1^s-1  exp  [-sx(£[l  +  S'^C'l'^~^s^b°d''ab]  +  sx(p/e)? 

(7) 
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(8) 


=  r1 


co 

i 

n=0 


1 

n! 


2 

m=0 


/  s(c+d-a-bj  + 
'  (s+a)(s+l 


cd- 
s+b) 


■ab\ni 


r(m-l/2) 

r(-i/2)r(m+i) 


Eq.  8  may  now  be  inverted  term  by  term.  It  should  be  noted 
that  each  term  in  this  equation  can  be  reduced  to  a  sum  of 
terms  of  the  form  of  a  constant  multiplied  by  s"1,  (s+a)_n, 
or  (s+b)"n  therefore  the  inverted  series  will  have  the  form 


F(t)  =  6(x)  +  e-atJQ  0n(x)t"  +  e~btJ0  !>n(x)tn  (9) 

Eq.  (9)  suggests  an  interpolation  of  the  form  of  Eq.  (1)  (with 
j  =  2)  for  any  given  value  of  x,  in  which  F(°o)  is  the 
value  of  the  response  as  t  ->  «  and  therefore  from  Eq.  (7) 

F(oo)  =  1.  The  sum  F(»)  -  A  -  A  represents  the  discontin- 

1  2 

uity  at  the  wave  front  and  may  also  be  obtained  from  Eq.  (7) 
using  the  usual  asymptotic  techniques  and  the  expansion  of 
(1+*)* 


F(t)  =  5<|;|^;grb  + 

I 

F(Q)  =  Lim  s  F(s)  =  exp[-x(p/e)2(c4d-a-b)/21 

s  ->  00 

=  F(oo)  -A  -A 
1  2 

Eq.  (1)  may  be  re-written 


,.]+sx(p/ 


(10) 


F  ( t )  ~1 


F(0)-1 
+  140 


+  e 


(ID 


where  a  =  A  /A  remains  to  be  determined. 
1  2 

of  Eq.  (8),  n  =  0,  is  s”1.  For  n  =  1, 


The  first  term 
there  is  the  term 
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(12) 


plus  terms  of  lower  order  in  s  as  s  ■+  ».  For  all  other 
n,  the  terms  of  Eq.  (8)  are  of  lower  order  in  s  as  s  ->  « 
than  Eq.  (12).  For  these  reasons  Eq.  (12)  is  used  to  deter¬ 
mine  the  constant  a.  The  inverse  transform  of  Eq.  (12)  is 


As  an  example,  let  G  =  G  =  1,  G  =  100,  x-1  =  C  =  1, 

2  3  1  1 

x"1  =  d  =  100.  It  follows  from  the  definition,  Eq.  (4),  that 
a  =  99.029,  b  =  .0099,  and  a  =  .961.  At  the  point 

i 

x(p/e)2x_1  =  1  on  the  rod,  the  interpolation  gives 


a  -  H(t-l) [l-.3l8(  .961  e-99-029(t-D  +  e-0.0099(t-l)n 

It  may  be  noted  that  quite  frequently  a  may  be  very 
much  larger  than  unity,  or  smaller  than  unity.  In  such  a 
case  only  one  of  the  time  constants  really  governs  the  re¬ 
sponse.  However,  it  does  not  seem  easy  to  predict  such  a 
situation  on  physical  grounds  prior  to  the  analysis. 
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